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Abstract 

Increasingly, stable isotope ratios of nitrogen (S^^N) and carbon (S^^C) are used to quantify trophic structure, though 
relatively few studies have tested accuracy of isotopic structural measures. For laboratory-raised and wild-collected plant- 
Invertebrate food chains spanning four trophic levels we estimated nitrogen range (NR) using 5^^N, and carbon range (CR) 
using 5^^C, which are used to quantify food chain length and breadth of trophic resources respectively. Across a range of 
known food chain lengths we examined how NR and CR changed within and between food chains. Our Isotopic estimates 
of structure are robust because they were calculated using resampling procedures that propagate variance in sample means 
through to quantified uncertainty in final estimates. To Identify origins of uncertainty in estimates of NR and CR, we 
additionally examined variation in discrimination (which Is change In 5^^N or S^^C from source to consumer) between 
trophic levels and among food chains. S^^N discrimination showed significant enrichment, while variation in enrichment 
was species and system specific, ranged broadly {1.4%o to 3.3%o), and Importantly, propagated variation to subsequent 
estimates of NR. However, NR proved robust to such variation and distinguished food chain length well, though some 
overlap between longer food chains Infers a need for awareness of such limitations. 5^^C discrimination was inconsistent; 
generally no change or small significant enrichment was observed. Consequently, estimates of CR changed little with 
Increasing food chain length, showing the potential utility of 5^^C as a tracer of energy pathways. This study serves as a 
robust test of Isotopic quantification of food chain structure, and given global estimates of aquatic food chains approximate 
four trophic levels while many food chains include invertebrates, our use of four trophic level plant-invertebrate food chains 
makes our findings relevant for a majority of ecological systems. 
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introduction 

Understanding food web structure is of critical importance to a 
broad suite of ecological theory given that trophic dynamics 
between individuals, populations, species and functional guilds 
underpin the ecological functioning and evolution of biological 
communities [1], [2]. Quantifying food web structure (trophic 
structure hereafter) is therefore a prerequisite to better understand 
how it in turn interacts with emergent properties of organisms and 
the environment, such as energy flux [3] , population dynamics [4] , 
patterns of biodiversity [5] and ecosystem functioning [6], [7]. 
Determination of feeding relationships between species is integral 
to quantifying trophic structure, and traditional methods include 
gut-content analysis, faecal analysis and behavioural observations. 
However, these methods can be laborious and may not reflect 
variation in digestibility and assimilation of source items, and if 
limited in their collection in space and time, may lead to over or 
under representation of source contributions [8]. Increasingly, 
stable isotope ratios of nitrogen (N''':N''', termed 5'''N) and 



carbon (C''':C'^, termed S'^C) in consumer tissues are used to 
provide a temporally and spatially integrated construct of dietary 
niche [8] , with 5 'N and 8 ' C of consumer proteins reflecting the 
proteins of their food sources [9], [10]. Typically, enrichment in 
5' 'N of 2.5%o to 3.4%o is observed from diet to consumer [1 1-13], 
allowing determination of an organism's trophic level [11], [14], 
[15] and overall food chain length [16], [17]. Conversely, 
enrichment in 5 C is much smaller between diet and consumer 
[11], [13], and because basal sources often dilTer in their S'^'C 
values, S'^^C can be used to trace prey-consumer connections or 
food chains [11]. Hence change in 5''^N and S'^C from source to 
consumer as described (termed trophic discrimination factors and 
represented as AS'''N or A8'''C), is the mechanism that crucially 
underpins the positioning of individuals, populations and species 
relative to one another in bivariate isotopic space (typically with 
5 'N on a y-axis and 5''^C on an x-axis). Importantly this 
subsequently allows for measures of Euchdean distances across the 
isotopic space occupied by populations, species or communities in 
order to quantify aspects of trophic structure [18-21]. For 
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instance, food chain length is calculated as X + (nitrogen range/ 
average; A8' 'N), where X is minimum trophic position and 
nitrogen rang(; (NR) is mean difference between trophic levels of 
maximum and minimum 5'''N [11], [14], [18]. Similarly, carbon 
range (CR) measures breadth of trophic sources and is calculated 
using mean difference between maximum and minimum 8''^C of 
community members [18]. Observational studies have largely used 
such measures to quantify food chain length, typically in aquatic 
systems, in response to factors such as ecosystem size, disturbance 
and productivity [22-24]. Thus the use of stable isotope ratios in 
an organism's tissues to provide temporally and spatially 
integrated dietary data is proving a very valuable methodology 
for trophic research. 

Critically though, variation in AS'^N and A6'"*C generates 
uncertainty/error in subsequent estimates of trophic structure and 
relationships [11], [17], [25-28]. Such variation in discrimination 
factors is well documented in the literature (— 0.8%o to 5.9%o for 
5'-'N and -2.7% to 3.4%o for 8''^C, excluding fluid feeders; [29]), 
This can be a consequence of multiple factors, including dietary 
protein quality, metabolic process and efficiency of protein 
assimilation and loss, fasting, growth rate, age, size, tissue type., 
sample size and sampling process, although there is considerable 
debate on which are most important [12], [13], [27], [30] (and 
references therein for all). It is the cumulative effect of all such 
underlying sources of variation that we observe as variation in 
discrimination factors between different individuals and species. 
Common practice often uses mean estimates of A5'^N and A5'^C 
ignoring variabilit)' around estimates; conserjuently, derived 
estimates of trophic structure and subsequent ecological conclu- 
sions may lack accuracy. 

These issues associated with variance have driven recent 
innovations in the analysis of isotopic data that provide 
practitioners with tools to apply Bayesian inference to the 
calculation of trophic structure metrics [20]. These approaches 
are ideally suited to testing effects of variance as they provide 
population and community trophic metric estimates based on 
resampling of variance in mean and 5'^C estimates, 

effectively quantifying and propagating variation in raw data as 
uncertainty in subsequent metric outputs, allowing for critical 
examination of precision in estimates of trophic structure. 

Despite variation in discrimination being widely documented 
(e.g. [11], [13], [29]) to date relatively few studies have validated 
isotopic structural measures such as NR and CR against known 
measures of food chain length, despite repeated calls for research 
[21], [27]. Of such studies, variation in baseline S'^'N has been 
shown to affect food chain length estimates [17], and variation in 
A6'^N to affect trophic level/food chain estimates [1 1], [25], [28], 
all in aquatic systems. Such studies have therefore been vital in 
guiding our understanding and use of these techniques; however, 
most have used wild systems which may not accurately identify 
basal sources or preclude omnivory, and thus do not explicidy 
know/control trophic levels against which to compare isotopic 
measures of structure. To further improve our understanding of 
isotopic approaches to quantifying trophic structure, so it would be 
prudent to additionally test the accuracy of isotopic measures of 
trophic structure for multi-trophic level food chains when trophic 
levels are explicitly controlled and food chain lengths known, in 
addition to allowing for the propagation of variation in A8' 'N and 
A8''^C to final isotopic trophic measures. 

It is of further importance to understand effects of variation in 
AS'^N and AS'^C upon structural measures such as NR and CR 
given these univariate measures underpin other bivariate (8''^N 
with 5''^C) measures of trophic structure in isotopic space. NR is 
also the most used isotopic metric in observational studies [21], 



and its function as a tool to quantify trophic level, food chain 
length or as a component of bivariate measures is dependent upon 
an assumed constant 5'^N enrichment with each consumer level. 
Variance around this assumed average enrichment constitutes 
unknown error in estimates of NR in observational studies. Thus, 
experimental validation would improve understanding of the 
importance of variance in A8'^N and A8'^C for affecting 
quantifications of trophic structure in wUd systems [21], strength- 
ening subsequently derived ecological conclusions, in addition to 
further catalysing the development and use of these techniques to/ 
b)' a ^vider audience of ecologists. 

In this study, we use natural plant and insect food chains raised 
under controlled conditions to examine the dynamics of 5'^N and 
S'^C over four trophic levels, and using resampling procedures to 
allow propagation of variance in AS'^N and AS'^C, test the 
accuracy of the isotopic metrics NR and CR to quantify trophic 
structure. We use a terrestrial plant-insect system in this study 
simply because terrestrial systems are currently less well studied 
relative to aquatic systems and it would therefore broaden our 
knowledge of isotopic techniques to have a greater range of 
validation studies to draw upon. By using four trophic levels, we 
broaden cross-system applicability of our results to a larger 
repository of aquatic studies, given that global aquatic food chain 
lengths have been estimated at 3.5 to 4.0 trophic levels [17]. 
Specifically, we test three questions: 1) How consistent are 
and A5'^C between trophic levels and among food chains? 2) Does 
NR accurately determine food chain length? 3) How does CR 
change with food chain length? 

Methods 

Ethics Statement 

AU wild animals used in this experiment are not protected by 
law, and were collected on land belonging to the University of 
Exeter for which permissions to collect had been granted. 

Experimental Set-up 

To test dynamics of 8''^N and 8'''C with changing trophic level, 
three replicate food chains were raised in the laboratory, with a 
further analogous food chain colkx tcd from the wild to allow 
comparison with wild systems. AU food chains had four trophic 
levels consisting of: primary producer (plant) — » herbivore (aphids 
feeding on plant phloem sap) predator (hoverffy larvae feeding 
on aphids) secondary predator (parasitoid hymenoptera which 
are obligate endoparasites of hoverffy pupa, with a single 
parasitoid emerging from a single pupa). These organisms were 
used as they provided us the opportunity to study natural food 
chains (that can be found in nature) and were easy to culture in 
controlled laboratory conditions so that trophic levels were 
explicitly known. Additionally, using natural food chains allowed 
us to directiy contrast observations of controlled lab raised food 
chains with our analogous wild collected food chain in order to 
interpret and ensure the relevance of our laboratory findings for 
drawing conclusions and lessons that could inform future studies 
using wild systems. 

For laboratory food chains, grain aphids [Sitohion avenai) were 
raised on two independent food plants; One based on a C3 
photosynthetic pathway (wheat Triticum ae.stivum) and the other 
based on a C4 pathway (maize ^ea mays), enabling separation of 
the plants on a S'^^C axis and thus broadening the generality of 
any observed patterns. Plants were raised on a common source of 
homogenised compost and distilled water, and introduced to 
aphids at 20 days (wheat) or 30 days (maize). T' generation larvae 
of wild-caught hoverfly [Syrphus vitripennis) were fed either wheat or 
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maize raised aphids or an approximate 50:50 ratio of both. Within 
each treatment, 24—48 hours after hatching, a random subset of 
hoverfly larvae were exposed to wild caught adult female 
parasitoids [Diplazon laetorius) to allow parasitic oviposition. All 
plants and insects were raised under a 16:8 light: dark cycle at 70% 
humidity. Plant leaves and aphids of all ages were collected at 
random and frozen (— 20°C) prior to tissue preparation. Hoverfly 
larvae entered pupation 8-10 days after hatching and after 
72 hours pupation were frozen (— 20°C) for later dissection. Prior 
experimentation identified 72 hour pupation as suitable to provide 
soft pupa tissue comparable to that likely consumed by parasitoid 
larvae (i.e. after exuviae had formed). Parasitised hoverfly larvae 
were allowed to complete pupae development (19-21 days) and 
newly eclosed adult parasitoids were frozen (— 20°C) within 
12 hours without having fed. 

For the wild food chain, nettle (Urtica dioica) leaves and nettle 
aphids (Microlophium carnosum) were collected independently and 
frozen (— 20°C) for later preparation. Hcnerfly larvae [Syrphus 
vitripemis) were collected when judged at >50% grown and 
laboratory raised on daily-collected wild nettle aphids until 
pupation, under a 16:8 light:dark cycle at 70% humidity. Pupation 
proceeded until either adult hoverflies or adult parasitoids 
(Diplazon laetorius) eclosed after 10-11 or 16-20 days respectively, 
and were frozen (— 20°C) within 1 2 hours without having fed. Our 
four replicate food chains are hereafter termed after their plants as 
wheat, maize, wheat + maize (w+m) and nettle. Parasitoid 
hymenoptera are referred to as wasp hereafter. 

Tissue Preparation & Lipid Extraction 

We used tissues for sources to best represent assimilation in 
consumers, as shown to be important [31]. While aphids fed on 
plant sap, we used whole leaf tissue given difflculties of extracting 
sap and because whole leaf tissue 5'''N has been shown not to 
differ from sap [32] . For each food chain, following dissection, soft 
internal tissues of 60-80 aphids were pooled to produce a single 
sample. Hoverfly larvae soft tissue was obtained from pupae 
casing. Wasps represented end consumers and we used whole 
tissues. For wild hoverflies and wasps we used adult whole tissues. 
Individual hoverflies and wasps each provided single rephcates. 
Sample sizes for plants, aphids and hoverflies were n= 15 (except 
nettie hoverflies, n = 7) while wasps were more difficult to obtain: 
wheat {n= 10), w+m (n = 7), netde (n = 6). For maize, no wasps 
were obtained due to high larval mortality. AH samples were dried 
at 45°C for >48hrs and homogenised. Subsequendy, insect 
samples were immersed in 2:1 chloroform: methanol solution for 
50 minutes to remove free lipid, and then left to air dry. Prior 
experimentation showed lipid extraction to have no significant 
effects on 5'*N of these insects. 

Stable Isotope Analysis 

For all samples, 0.5 mg±0.05 (insect) or 3 mg±0.1 (plant) dried 
material was enclosed in tin [:apsules. Stabk- isotope analysis (SLA) 
was conducted at the Food and Environment Research Agency, 
York, UK. Samples were analysed for and in a Fisons 
EAl 108 elemental analyser (Carlo Erba Instruments, Milan, Italy), 
coupled with an Isoprime isotope ratio mass spectrometer (GV 
Instruments, Manchester, UK). Stable isotope ratios are reported 
in delta (8) notation where 5''''N and 8''^C = [(Rsam],i< /Rsiiin<iard) 
-1] xlOOO, where R is '-'N/'^N or '''C/'^C. Isotope ratios are 
expressed in per mil (%o) relative to the ratio of international 
reference standards (Rstandani) which are Atmospheric Nitrogen 
and Vienna PeeDee Belemnite (VPDB) for nitrogen and carbon 
respectively. Measures of standards placed throughout samples 
exhibited acceptable instrument reproducibility of <0.09%o (SD) 



for and <0.18%o (SD) for S''^C using collagen standard, 

insect whole tissue standard (cockroach; Nauphoeta cinerea), and 
sucrose C4 plant standard. 

Data Analysis 

Initial analyses used two-way analysis of variance (ANOVA) to 

test effects of the two explanatory variables trophic level (levels = 
plant, aphid, hoverfly or wasp) and food chain type (levels = wheat, 
maize, w+m or nettle) on 8''^N and then 8'''C. Model 
simplification used backwards stepwise regression from a maximal 
model and ANOVA model comparisons to identify non-significant 
terms for elimination. Homogeneity of variances and normality of 
model residuals were checked in all instances. To determine where 
significant differences lay between levels within treatments, 
subsequent one-way ANOVA for each food chain type were 
tested with Tukey post hoc tests, for both 8'''N and S''^C. 

Within each food chain type, we then calculated mean (± SD) 
AS'^N and AS'^C as the difference in 8''^N or 8''^C between each 
source and its consumer by randomly pairing replicates (n = 6 to 
15). To establish underlying sources of variation in AS'^N, two- 
way ANOVA tested effects of explanatory* variables trophic link 
(levels = plant-aphid, aphid-hoverfly, or hoverfly-wasp) and food 
chain ype (levels = wheat and nettle). Only nettle and wheat food 
chain data were used as these contained all three trophic links, 
allowing for a balanced analysis. For explanatory variables, 
variation in AS'^N was quantified using sums of squares in model 
outputs and was expressed as a proportion of the nuU model 
variance. Model simplification was used as described above. AU 
analyses were conducted in R version 2.14.1 [33]. 

Nitrogen range (NR) and carbon range (CR), calculated as the 
mean diflerence between trophic levels of maximum and 
minimum 5'^N or respectively, are quantifications of trophic 
structure with NR representing food chain length and CR the 
breadth of energy sources. For all species combinations of each 
food chain length within each food chain type separately, NR and 
CR were independentiy calculated using resampling of uncertainty 
around sample mean estimates to provide probabilistic distribu- 
tions representing 50%, 75% and 95% credible intervals of mean 
estimates for NR and CR, using the SIBER computational code 
[20] in the R package SIAR [34], [35]. Additionally for NR, for 
each food chain length, resampled mean estimates across food 
chain types and species combinations were pooled to produce 
overall 95% credible intervals. 

As established by preliminary experiments, we applied correc- 
tion factors of-O.7%0 for 8'"^N and +0.4%o for 8'''C to hoverflies 
on the wild nettie food chain (which used adult tissues) to make 
them directiy comparable to larval hoverflies on laboratory food 
chains. 

Data Availability 

All raw data used in this study are provided as supplementary 
material (Dataset SI). 

Results 

How consistent are AS^^N and A8^^C between trophic 
levels and among food chains? 

8^^N. Though given the same nitrogen source, 8''N of wheat 
plants were depleted relative to maize plants by ~2%o, while wild 
nettie plants (of independent nitrogen source) were slightiy 
enriched (0.3%o) relative to maize. Such difiTerences between food 
chain types were largely propagated to higher trophic levels 
(Fig. la) and found to be significant; two-way ANOVA showed a 
significant interaction between trophic level and food chain type 
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Figure 1. S N (a) and S C (b) across four trophic levels of four replicate terrestrial food chains. Mean ± SD (%o) are shown. n = 6 to 15. 
For plant S^^N, mean ± SD are offset on x-axis for clarity. Dotted lines are trophic links between two aphid prey sources and their hoverfly predator. 
doi:1 0.1 371 /journal.pone.0093281 .gOOl 



indicating efTects of these variables were interdependent (Fjg 152) 
= 11.57, P<0.001). To determine if large variances in 8''^N of 
maize and wheat plants (Fig. la) disproportionately affected these 
results, we repeated this analysis excluding these plants, but found 
no dilference in outcome. 

To ascertain patterns of 5''^N discrimination, Tukey post hoc 
tests performed on one-way ANOVA for each food chain type 
estabhshed where 5''^N differed between trophic levels. Overall, 
discrimination factors showed significant enrichment from 
source to consumer (range of 1.4%o to 3.3%o), in all but three 
instances (Table 1). Exceptional to this trend, wheat feeding aphids 
showed significant average depletion in 5''^N (— 2.4%o) relative to 
hosts. 

Given enrichment was broadly consistent and larger than 
that observed in 5''^C (Table 1), for 5''^N we also used two-way 
ANOVA to determine sources of variation in AS'^N (Table 2). 
After disproportionate plant variation was excluded, we observed a 
significant interaction between trophic link and food chain type 



accounting for 32% of variation in A5''N. This significant 
interaction infers variation in AS'^'N as caused by these variables 
was species and system specific. 

6^^C. 5'''C values were ~17%o different between wheat and 
maize food chains, with hoverfly and wasps on the w+m food 
chain approximately half way between the two having integrated 
aphid sources from both (Fig. lb). Two-way ANOVA showed 
5''^C was significantiy affected by an interaction between trophic 
level and food chain type (Fj^ 152) = 19.22, P<0.001), indicating 
food chain effects on 5''^C were affected inconsistently by trophic 
level (Fig. lb). 

Differences in 5''^C discrimination across trophic links and food 
chain types were of variable direction and magnitude (— 0.7%o to 
1.9%o, excluding w+m aphid-hoverfly; Table 1). Across trophic 
links and food chain types, Tukey post hoc analyses showed either 
significant enrichment (0.6%o to 1.9%o) or no change in 5 ' C 
between trophic levels (Table 1). 



Table 1. Mean ± SD (%o) A5^^N and AS^^C across four trophic levels of four terrestrial food chains. 





Source 




Consumer 


Food Chain 














wheat 


maize 


w+m 


nettle (wild) 


Plant 




Aphid 


-2.40±2.63«** 


1.12±2.58 




1.40±0.75*** 


Aphid 




Hoverfly 


2.1 8± 1.25** 


1.51 ±0.69* 


4.81 ±0.87***(w) - 
0.92±1.02***(m) 1.95±3.06* (c) 


3.08±0.60*** 


Hoverfly 




Wasp 


3.31 ±1.00*** 




2.49±0.62*** 


1.38±0.61** 


S"C 


Plant 




Aphid 


0.65 ±0.92** 


1.94±0.36*** 




1.88±0.86*** 


Aphid 




Hoverfly 


-0.04±0.33 


-0.69±0.42*** 


9.33±0.87***{w) - 
9.05±0.79***{m) 0.14±9.38 (c) 


1.22±0.71*** 


Hoverfly 




Wasp 


0.20±0.30 




0.39±0.65 


0.40±0.31 



Significant differences in 6^^N and 5^^C between source and consumer within each food chain indicated as P= <0.05*, <0.01**, <0.001***. w+m hoverfly has two aphid 
sources and values are given for both: (w) = wheat aphids, (m) = maize aphids, (c) = combined. 
doi:1 0.1 371/journal.pone.0093281 .tOOl 
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Table 2. Sources of variance in AS^^N (%o). 







All combinations 


Excludes primary producers 


Null model variance in A8'^N 


398.1 


52.8 


Model Term 


Proportion of null variance explained 


Proportion of null variance explained 


Food Chain 


0.04 


0.02 


Trophic Link 


0.42 


<0.01 


Food Chain * Trophic Link 


0.21 **• 


032 


Full model (all terms combined) 


0.67 


0.34 



Variance in AS'^N as accounted for by either food chain or tropic link was established by expressing ANOVA model terms as a proportion of the null model variance. 
A6^^N is based upon differences in raw 5^^N between trophic levels, for all source-consumer links on wheat and nettle food chains only (to provide a balanced analysis). 
ANOVA was conducted twice: firstly for all source-consumer links and secondly excluding all links including a primary producer. A significant interaction explained 
greater variance in AS^^N after primary producer links were excluded. Significant model terms are indicated as *** with P = <0.001. 
doi:l 0.1 371 /journal.pone.0093281 .t002 



Does nitrogen range accurately determine food chain 
length? 

In most instances, estimates of nitrogen range (NR) were 
observed to increase with greater food chain length within each 
food chain (Fig. 2). Exceptional to this were low NR estimates on 
the wheat and w-l-m food chains for combinations including wheat 
plants. 

Within a food chain, in some instances modal or credible 
interval values of NR varied distinctly between different species 
combinations of the same food chain length, but such differences 
were not constant across food chains suggesting a species and 
system specific nature to such variation. Across food chains, 95% 
credible interval estimates of mean NR were not generally larger 
than l%o 8''^N either side of the modal value (Fig. 2). The 
exception was larger estimates on the wheat food chain when 
wheat plants were included as a consequence of uncertainty in 
mean estimates (Fig. 2), propagated from large sample variation in 
wheat plants (Fig. 1). 

Overall estimates of NR based on combining resampled mean 
estimates from all four replicate food chains and all species 
combinations (for each level of food chain length), showed modal 
NR values to increase by between 1.2%o and 2.7%o with each 
additional trophic level (Fig. 3). Overlap in 95% credible intervals 
between different food chain lengths (Fig. 3a) was reduced when 
wheat plant combinations were excluded (Fig. 3b), and then 
further reduced when maize plant combinations were also 
excluded (Fig. 3c). Notably, there was no subsequent change in 
95% credible intervals when nettle plants were additionally 
removed (Fig. 3d), suggesting that these wild plants did not 
contribute noticeably to variation in estimates of NR. Following 
removal of wheat and maize plant combinations, increase in 
modal NR with trophic level was more closely matched at 1.7%o 
and 2.4%o. Removal of nettle plant combinations increased modal 
values, most likely due to nettle-aphid discrimination factors being 
lower than some other trophic links (Table 1). 

How does carbon range change with food chain length? 

For estimates of carbon range (CR), we observed few consistent 
patterns in CR across differing food chain lengths within or among 
food chains (Fig. 4). Wheat and maize food chains both showed no 
pattern of change in CR with changing food chain length; wheat 
CR was <l%o across two, three and four trophic levels, while 
maize CR was 2%o for combinations of both two and three trophic 
levels. For the w-l-m food chain there was marginal increase in 
modal CR with food chain length (two to four trophic levels 
~0.3%o to 2.3%o) but this was considerably less than variation in 



mean estimates as shown by large 95% credible intervals (>3%o), 
on account of two isotopicaUy disparate 5''^C plant sources. 
Conversely, the wild netde food chain trended towards greater 
modal CR with food chain length (two to four trophic levels ~ l%o 
to 3.5%o); however, overlap in 95% credible intervals between 
different food chain lengths was observed, while modal values of 
different species combinations of the same food chain length also 
differed by =l%o to 1.5%o. 

Discussion 

We tested the accuracy of isotopic measures of trophic structure 
against known trophic positions using well replicated terrestrial 
plant and invertebrate food chains spanning four trophic levels. 
We found that despite some overlap in nitrogen range (NR) 
between longer food chains, across a range of different species 
combinations and food chains NR generally quantified food chain 
length well, suggesting robustness to observed variation in 
discrimination. Additionally, we found few consistent trends in 
5''^C discrimination with typically small (0.6%o to 1.9%o) or no 
enrichment, and concurrendy litde and inconsistent change in CR 
with food chain length, emphasising the utility of 5'^C to trace 
diet-consumer pathways. We suggest our estimates of food chain 
trophic structure are particularly robust because they were 
calculated using resampling procedures, allowing for propagation 
of variance in AS'^'N and AS'^'C into quantified uncertainty in 
final structural estimates. 

How consistent are A§^^N and AS^^C between trophic 
levels and among food chains? 

6*^N. Within species variation in 5'-'N of 0.3%o to 0.9%o 
(excluding plants) was observed, and was a likely consequence of 
differences in AS'^N between individuals due to compositional 
differences in food consumed and metabolic differences in 
assimilation of that food [13], [27], [36] (and references therein). 
An additional anomalous observation was large 5'^N variation in 
wheat and maize plants that was not observed in wild nettle plants, 
suggesting that laboratory conditions affected variation. As such 
variation was not observed for 5"C which is sourced from the 
atmosphere and showed much smaller comparable variation 
across wheat, maize and netde plants, we speculate 5''^N variation 
was explained by micronutrient pockets in laboratory soil medium 
through incomplete homogenisation. Alternately, stress, as caused 
by unnaturally high and variable laboratory aphid density, may 
have affected plant metabolism and hence nitrogen balance, as 
shown for other taxa [37], [38]. As we pooled aphids to produce 
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Figure 2. Probability distributions of mean nitrogen range (NR) for different food chain lengths for each food chain. NR Is difference 
between mean S^^N (%o) of community end members, based on resampling (n = 10,000) of uncertainty In sample mean estimates. Black dots 
represent mode (of means), while shaded boxes (dark to light) show 50%, 75% and 95% credible Intervals for mean estimates. No wasps were 
obtained for the maize food chain. x-axIs labels are species identity: p = plant, a = aphid, h = hoverfly, w = wasp. Parenthesis number shows food chain 
length. 
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each aphid sample, plant 5' 'N variation likely averaged across 
these multiple aphids, explaining why large 5 'N variation was not 
subsequently seen in aphids. 

We found significant 5''^N enrichment from source to consumer 
in a majority of instances, varying between 1.4%o and 3.3%o, a 
range concurrent with literature estimates [11-13]. Averaged 
across food chains, such enrichment was 2.2%o (excluding 
discrimination from wheat and maize plants to their aphids). This 
value is marginally lower than literature average estimates for 
invertebrates of 2.5%o [12], [13], and lower than the overall 
literature estimate of 3.4%o which is commonly employed to 
calculate trophic levels and food chains (e.g. [1 1], [23], [24]). This 
higher overall literature value may be reflective of trophic 
discrimination factors associated with vertebrates [13] and thus 
the lower enrichments noted in our insects re-emphasises the need 
to use taxa-specific discrimination factors when using isotopic data 
to calculate ecological parameters [13]. Additionally, broad 
enrichment range of 1.9%o (1.4%o to 3.3%o) in our results, as 
concurrent with previous studies [11-13], demonstrates variation 
in trophic discrimination factors is an important artefact in 
isotopic data that should be accounted for. Thus the use of 
averaging to produce commonly used trophic discrimination 
factors excludes variation in discrimination from final ecological 



estimates. Given such variation was present in our system, our 
subsequent use of resampling procedures allowed us to propagate 
discrimination uncertainty into our estimates of food chain 
structure, producing estimates that we contend are more accurate 
and hence ecologically robust [20]. 

In this study, as a single exception in source to consumer 
enrichment, wheat feeding aphids were on average depleted in 
5' 'N relative to hosts, concurring with other studies [32], [39]. A 
previously observed negative relationship between aphid 
discrimination and host plant total nitrogen content [32] is 
consistent with our results; aphid enrichment was observed on 
maize and nettie plants that had low total nitrogen contents of 
1.7% and 2.7% respectively, while wheat plants had higher 
nitrogen of 5.4%. 

While examining sources of variation in A5' 'N, we found a 
significant interaction between trophic link and food chain type, 
suggesting effects of these variables were inconsistent and thus 
could not be generalised. This result is in close agreement with 
recent observations of discrimination variation at broad spatial 
scales in wild systems whereby variation was dissimilar across taxa 
within and between trophic groupings [28]. Mechanistic under- 
pinnings explaining variation in A5''^N are beyond the scope of 
this study but briefly, for different trophic links and for the same 
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Figure 3. Overall nitrogen range (NR) for food chains containing two, three or four trophic levels. Values derived from combining 
resampled mean values from all four replicate food chains. TL = trophic levels. Circles represent mode (of means) and bars 95% credible Intervals for 
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trophic links on difiFerent food chains, variation in A5' '^N was likely 
due to differences in nitrogen [32] and protein [40] quality of food, 
and/ or consumer's differing metabolic assimilation of that food 
[27] (and references therein). 

6^^C. Trophic discrimination of 5 ' C from source to con- 
sumer was inconsistent between trophic links and food chain type, 
ranging from — 0.7%o to 1.9%o, and more generally observed as 
either showing no change or significant enrichment of between 
0.6%o to 1.9%o, concurring with discrimination reviews [1 1], [13]. 
Given that in our system the carbon axis was broad, such 
inconsistent and generally small trophic discriminations meant 
was diagnostic of food chain type. Similarly, S'^^C signatures 
of hoverfly predators on the w+m food chain were intermediate of 
their two disparate aphid sources, illustrating well the usefulness of 
5''^C data to integrate and reflect dietary sources. 



Does nitrogen range accurately determine food chain 
length? 

Using resampling procedures to calculate nitrogen range (NR), 
we tested how NR changed with known food chain length within 
food chains, and then across food chains types. Excluding wheat 
plant-aphid combinations, we found that NR accurately deter- 
mined food chain length within all food chains. This finding is 
concurrent with a previous validation showing positive correlations 
between trophic positions of freshwater fish estimated by both 
5'''N and traditional gut content analysis [14]. 

In our study, the inclusion of wheat plant-aphid combinations 
depressed NR measures on the wheat and w+m food chain types. 
This was because aphids were depleted relative to wheat plants, 
such that aphids were effectively base of the S''^N food chain while 
wheat constituted an additional trophic level that did not act to 
extend NR. This importantly shows exceptional species-specific 
effects may adversely affect the accuracy of isotopic measures of 
food chain length. As previously we had identified the wheat 
plant-aphid relationship as an exception to the generic enrich- 
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Figure 4. Probability distributions of mean carbon range (CR) for different food chain lengths within each food chain. Distributions 
based on resampling (n = 10,000) of uncertainty in sample mean estimates. 50%, 75% and 95% credible intervals as in Fig. 2. x-axis: p = plant, 
a = aphid, h = hoverfly, w = wasp. Numbers in parenthesis show food chain length. 
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ment in 5' 'N from source to consumer, we feel justified in 
concluding that, more generally, NR predicted food chain length. 

Overall, when all trophic combinations for all food chain types 
were combined, modal NR values increased by between 1.2%o and 
2.7%o with each trophic level (or 1.7%o and 2.4%o excluding wheat 
and maize plant combinations), suggesting a robustness of this 
technique for calculating food chain length. Excluding combina- 
tions including wheat and maize plants (on account of large 
variation in their 5 N signatures), 95% credible intervals of NR 
estimates showed some overlap between food chain lengths of 
three and four trophic levels, inhibiting accurate estimation of food 
chain length at these points of overlap. Similarly, for wild-sampled 
freshwater algae, invertebrate and fish food chains it has been 
shown that differences in 5'^N between trophic levels was largely 
reflective of food chain length, though with variation in A5''^N 
causing some uncertainty/ error in estimating trophic level [28]. 
Current practice uses NR to estimate food chain length and 
determine subsequent conclusions, but rarely has NR's use been 
tested, with terrestrial systems particularly understudied [27]. 
Based on our results, we urge caution interpreting food chain 
length when NR values fall in known overlap boundaries; i.e. in 
our study NR values of 5%o or 6%o could be either three or four 
trophic levels. Such overlap was a consequence of variation in NR 
values within each trophic link and between trophic links of the 
same food chain length, both within and between food chain types. 



as ultimately caused by variation in A5'^N within and between 
source-consumer pairings (Table 1). Of the few studies that have 
examined effects of variation in baseline 5' 'N [17] and AS''^N 
[11], [25], [28] on error of trophic level or food chain length 
estimates, all have found error concurrent to the variation we show 
in estimates of NR. Importantly, our research diversifies these 
studies by testing empirical rather than mainly theoretical 
measures of food chains, in a controlled laboratory context, using 
understudied terrestrial and invertebrate systems, in addition to 
providing robust estimates of NR based on resampling of 
variation. 

Additionally, as different food chain types in our experiment 
differed in their NR values for given food chain lengths as a 
consequence of variation in A5''^N, so it is worth urging some 
caution when using NR for direct comparisons between systems of 
simple single-pathway food chains. It has been reasonably 
suggested that such variation may average out over multiple 
trophic levels or larger sample sizes [27], though based on our 
results we would call for research to further test the importance of 
such variation. Pragmatically, as in our study the effect of food 
chain type on A5''^N was variable, then studies sampling across 
multiple food chains will absorb such variation as 'noise' that may 
be acceptable or specifically filtered out if food chain identity was 
categorised in analyses as a random effect. We speculate that such 
variation may also be averaged out when considering overall 
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structure of larger food webs containing multiple food chain 
pathways. 

Our overall estimates of NR were based on exact trophic 
positions in well replicated food chains and used resampling 
procedures that propagated uncertainty in discrimination factors 
to produce estimates that represented a full range of possible NR 
values. Given such an approach is likely to produce robust 
estimates of NR, and that more generally these estimates suggested 
that NR can accurately discriminate different food chain lengths, 
we conclude NR is a useful isotopic metric for quantifying food 
chain length. 

How does carbon range change with food chain length? 

Using resampling procedures we calculated carbon range (CR) 
values for all food chain length combinations on all four food chain 

types. Overall, our results show that across trophic levels CR 
changed litde and inconsistently both within and between food 
chain types (Fig. 4). Little change in CR over four trophic levels 
(modal range across food chain types <l%o to 3.5%o; Fig.4) 
suggests fidelity of values between primary producers and top 
predators, as concurrent with our earlier findings of small and 
inconsistent 5 C discrimination (Table 1). Given global estimates 
of aquatic food chains approximate 3.5 to 4.0 trophic levels [17] 
htde change in CR over four trophic levels in our study 
demonstrates more broadly the utility of 5'^C to trace energy 
pathways, be they either simple singular chains or potentially when 
embedded within larger trophic structures. However, a potential 
caveat of using 5''^C is that source signatures may need to be 
distinct in order to determine different energy pathways; we 
acknowledge that our S'^C axis was broad (~17%o) with diflFerent 
plant species showing separation in S'^'C values. In systems where 
multiple consumers feed on a single basal source, or basal sources 
have closely aligned S'^'C, identification of different energy 
pathways may not be possible. 

Stable isotope approaches to quantifying food web structure 
continue to proliferate in ecological studies and so improving our 
understanding of the accuracy and limitations of these techniques 
is of importance. In this study, source-consumer 5'^N discrimina- 
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